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Dilute or semi-dilute solutions of non-intersecting self-avoiding walk (SAW) polymer chains are 
mapped onto a fluid of "soft" particles interacting via an effective pair potential between their centers 
of mass. This mapping is achieved by inverting the pair distribution function of the centers of mass 
of the original polymer chains, using integral equation techniques from the theory of simple fluids. 
The resulting effective pair potential is flnite at all distances, has a range of the order of the radius 
of gyration, and turns out to be only moderately concentration-dependent. The dependence of the 
effective potential on polymer length is analyzed in an effort to extract the scaling limit. The effective 
potential is used to derive the osmotic equation of state, which is compared to simulation data for the 
full SAW segment model, and to the predictions of renormalization group calculations. A similar 
inversion procedure is used to derive an effective wall-polymer potential from the center of mass 
density profiles near the wall, obtained from simulations of the full polymer segment model. The 
resulting wall-polymer potential turns out to depend strongly on bulk polymer concentration when 
polymer-polymer correlations are taken into account, leading to a considerable enhancement of the 
effective repulsion with increasing concentration. The effective polymer-polymer and wall-polymer 
potentials are combined to calculate the depletion interaction induced by SAW polymers between 
two walls. The calculated depletion interaction agrees well with the "exact" results from much more 
computer-intensive direct simulation of the full polymer- segment model, and clearly illustrates the 
inadequacy - in the semi-dilute regime - of the standard Asakura-Oosawa approximation based on 
the assumption of non-interacting polymer coils. 
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I. INTRODUCTION 



Polymer solutions have attracted the attention of 
theorists and experimentalists alike for many decades, 
and a theoretical understanding of their structure and 
phase behavior, based on mean-field and on scaling or 
renormalization group arguments, is by now well estab- 
lished 0-0| . Recently, there has been a growing interest 
in the structure, phase behavior and rheology of binary 
systems involving colloidal particles and non-adsorbing 
polymer |7|-|2^]. In such mixtures the mean size of the 
polymer coils, i.e. their radius of gyration Rg, is compa- 
rable to, or smaller than the diameter a of the colloidal 
particles. Since the latter may, for most purposes, be 
modeled as "hard" convex bodies dominated by excluded 
volume effects on the mesoscopic scale cr, it is clear that a 
statistical description of the polymer coils requires a high 
degree of coarse-graining to provide a tractable theory of 
these mixtures. Such coarse-graining is, more generally, 
desirable for theoretical investigations of large scale phe- 
nomena involving large numbers of interacting polymer 
chains in the dilute of semi-dilute regimes. In particu- 
lar, simulations of solutions involving many interacting 
polymer chains become rapidly intractable if a detailed 
description at the level of monomers or even of Kuhn 



segments is retained. It is therefore tempting to con- 
sider polymer coils as "soft" particles, and to replace the 
detailed interactions between segments by an effective 
interaction acting between the centers of mass (CM) of 
different polymer coils as shown schematically in Fig. |l|. 
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FIG. 1. Modelling polymer coils by effective "soft par- 
ticles". The A*' polymers, each made up of L segments, are 
replaced by A'^ particles interacting with an effective pair po- 
tential. The centers of the particles correspond to the poly- 
mer CM. The interaction of the polymers with a hard wall is 
modelled by a single soft particle-wall interaction. 



Similarly, an effective interaction must be worked out 
between the "soft" polymer coils and the "hard" colloidal 
particles. Such a drastic reduction in the number of de- 
grees of freedom, achieved by formally averaging over the 



coordinates of individual polymer segments, leads to a 
considerable simplification of the initial problem involv- 
ing Nc colloidal particles and NpL polymer segments, 
where Np is the number of polymer coils and L the num- 
ber of monomers or segments per polymer (i.e. the length 
of a polymer). The idea of representing a polymer coil 
by a single particle of radius of the order of Rg goes back 
to the work of Flory and Krigbaum pS], who consid- 
ered the infinite dilution limit of two isolated interact- 
ing polymers. A brief outline of subsequent theoretical 
and numerical work on the two-coil problem is given in 
Section [II. In this paper we generalize the idea to fi- 



nite concentrations, i.e. to dilute and semi-dilute polymer 
solutions. The effective interaction between the CM of 
polymer coils is determined by a combination of Monte 
Carlo (MC) simulations of a detailed segment model of 
interacting polymers, and of an inversion technique which 
allows the effective pair interaction to be extracted from 
the MC results for the center of mass pair distribution 
function. A similar inversion technique is applied to the 
density profiles of the CM of the polymers near a hard 
wall to determine the effective interaction between a wall, 
impenetrable to the polymer segments, and the CM of in- 
teracting polymers. The effective polymer-polymer and 
wall-polymer interactions provide a first step towards a 
complete description of colloid-polymer mixtures, with 
the hard wall considered in this paper representing a sin- 
gle colloidal particle of infinite radius. The ultimate goal 
is to go well beyond the familiar Asakura-Oosawa (AO) 
model which considers polymers to be non-interacting 
point particles, excluded from a sphere of radius a/2 + Rg 
around each colloidal particle |23|. This model leads to 
the well-known AO depletion interaction between hard 
sphere colloids [p4 - 26 1 . As an application of the general 
method outlined in this paper, the limitations of the AO 
picture will be illustrated in a calculation of the depletion 
interaction between two parallel hard walls. The effective 
interaction between polymer coils will be shown to lead 
to considerable deviations from the AO results, even in 
the dilute regime. 

A preliminary account of parts of the present work 
has been published elsewhere p^. A related soft parti- 
cle picture has recently been applied to polymer melts 
and polymer blends |28[ . However, the phenomeno- 
logical coarse-graining procedure proposed by these au- 
thors, and its practical implementation, differ consider- 
ably from the present "first principles" approach, which 
is better adapted to dilute and semi-dilute polymer solu- 
tions. Both methods are good examples of current efforts 
to bridge widely different length and time scales in com- 
plex fluids. 



II. SIMULATION MODELS AND METHODS 



Many physical properties of polymers in solution al- 
ready emerge from simple models which ignore chemical 
detail and describe the polymers as self avoiding walks 
(SAW) with hard segments interacting through a simple 
potential. For example, solutions of linear polymers in 
a good solvent are well modeled by N athermal SAW's, 
each made up of L non-intersecting segments, on a cu- 
bic lattice of M sites, with periodic boundary conditions. 
This model captures the leading scaling behavior and has 
been used for many decades to describe polymer solutions 
||l|-^. Slightly more sophisticated models exist, such as 
the fluctuating bond model Hg] or off-lattice hard sphere 
chains |30|, but the SAW lattice model is simple, efficient 
and allows for comparisons with previous studies. 

Within the lattice model, the monomer packing frac- 
tion is equal to the fraction of lattice sites occupied by 
polymer segments, c = N x L/M, while the concentra- 
tion of polymer chains is pb ~ c/L = N/M. For a single 
SAW chain, the radius of gyration scales as Rg ~ L'^, 
where v ~ 0.6 is the Flory exponent |1|. The overlap 
concentration p* , signalling the onset of the semi-dilute 
regime, is such that ATTp*Rg/3 ~ 1, and hence p* ^ L~^^ 

13- 

To sample the configuration space of the polymer sys- 
tem we employ the Monte Carlo pivot algorithm [BO 32 
which attempts to rotate part of the polymer around a 
random segment (the pivot). If the new trial configura- 
tion shows no overlap, the move is accepted, otherwise 
the old configuration is restored. This simple scheme 
turns out to be very effective for single polymers and di- 
lute polymer solutions where we found that it efRciently 
samples configurational space up to densities Pb/ P* ~ 1 
for L = 500 polymers. Because the polymers are re- 
stricted to a cubic lattice, the pivot move can only take 
place in 5 possible directions. For efhciency we store the 
complete lattice in memory, so that overlap between dif- 
ferent polymers can be easily checked for. In this way 
one has only to check of order L sites per polymer move, 
which is much more efficient than the NL'^ sites needed 
when each pair of segments has to be tested for overlap. 

In addition to the pivot moves, we also attempt to 
translate the polymer. This Monte Carlo move enhances 
the relaxation to equilibrium of the polymer solution, 
although the acceptance ratio for this move decreases 
rapidly if the density exceeds Pb/p* ~ 1 (for L = 500 
polymers). For densities deep in the semi-dilute regime, 
Pb/ P* > 1, we therefore also perform configurational bias 
Monte Carlo (CBMC) moves [|||4), in which part of 
the interior polymer is regrown. In addition, we attempt 
reptation moves where a limited number of segments at 
one end of the polymer are removed and regrown at the 
other end. By regrowing the polymer a bias is introduced, 
which is then corrected for in the sampling Pq,pjl . In the 



simulations at high densities, we find that we can regrow 
groups of up to about 20-40 segments in a CBMC move 
with a reasonable acceptance ratio (about 40 - 50 %). 
More sophisticated algorithms for very dense polymer 
systems are available pS], but are not necessary in our 
relatively dilute systems. 



III. EFFECTIVE POTENTIALS: TWO ISOLATED 
POLYMERS 

The theory of the effective interaction between two 
polymer coils in dilute solution has a long history. The 
first calculations were by Flory and Krigbaum in 1950 
p3[ , who showed that, within a mean-field picture, SAW 
polymers in a good solvent have a strongly repulsive in- 
teraction of the form: 
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where r is the distance between the CM of the two poly- 
mer coils, X is the usual Flory parameter, and /? — 1/k^T 
is the reciprocal temperature, with fee denoting Boltz- 
mann's constant. As long as the polymers are in a good 
solvent, the chains can be regarded as athermal, and for 
that reason we set /3 = 1 in the rest of this paper. 

The interaction strength at full overlap, Wj (r = 0), 
can be understood from the following argument: each 
polymer coil has a density of monomers c ~ L/V ^ while 
the volume of a polymer scales as V (x L^"^, so that 
c ex L^^^'^ (here v w 0.6 is the Flory exponent). If two 
polymers overlap completely then the mean-field free- 
energy of interaction would be proportional to the num- 
ber of monomers times the probability of contact of two 
monomers on different chains: 
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which implies that the polymer repulsion increases with 
polymerization L, and is typically much larger than fc^T. 
In an elegant paper, Grosberg, Khalatur, and 
Khokhlov iQ showed that Flory's argument was in fact 
incorrect. From scaling theory it follows that the proba- 
bility of an interaction between two monomers on differ- 
ent chains scales as c^/^^'^~^^ ~ c^"^ |37| instead of simply 



so that the free energy of interaction scales as: 
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In other words, the free energy of interaction at full over- 
lap of two equal-length polymers is independent of the 
degree of polymerization; polymer coils are not nearly as 
"hard" as one might naively expect. 

Kriiger, Schafer, and Baumgartner pq] put these ideas 
on a firmer footing using elaborate renormalization group 
(RG) calculations. In particular they calculated the full 



free-energy of overlap of polymers as a function of the 
CM distance. They found V2{r = 0) = 1.53e from an 
r-space e expansion and V2{r = 0) = 0.94e + 0.62e^ from 
a fc-space e expansion (e = 4 — d, so e = 1 for 3 dimen- 
sions). Although there are still significant quantitative 
differences between an 0(e) and an C(e^) calculation, 
implying that the e expansion has not quite converged, 
the qualitative picture is clearly that of a repulsive Gaus- 
sian type potential, as shown in Fig. 0. These calcula- 
tions were confirmed by a number of computer simulation 
studies, notably those by Olaj and collaborators p9|, and 
by Dauntenhahn and Hall pOl . 




FIG. 2. Effective CM-CM pair potential V2{r) for two iso- 
lated SAW polymers, here shown for different lengths L. The 
X-axis is scaled with Rg, to allow comparison. The pair in- 
teraction V2 (r) is approximately Gaussian. The height of the 
potential at r = decreases with length. Also shown is the 
RG result from an order 0{e^) expansion. Inset: 'iTvr^V2{r), 
which is more relevant to the thermodynamics of polymer so- 
lutions, shows much less variation with length L than V2{r). 



We repeated the calculation of the effective interac- 
tion between two isolated SAW polymer coils, to make 
sure that the simulations are carried out under conditions 
sufhciently close to the scaling limit. For two polymers 
at infinite dilution, the effective interaction can be de- 
termined by calculating the normalized probability P{r) 
of finding their respective CM's at a separation r. The 
effective potential V2(r) is then defined as 



V2ir) ==-ln(P(r)) 



(4) 



In the course of the simulation wc sample configurations 
of two polymers infinitely far apart using only the pivot 
algorithm. After every 1000 pivot moves, we calculate 
the overlap probability as a function of CM distance, by 
moving the polymers towards each other while check- 
ing for overlap. In addition, the radius of gyration is 
calculated for each length considered, from L = 100 to 



L — 8000. This reproduces the well known Flory seal- Since the potentials scale with r/Rg, this means that 



ing law Rg ~ L" . The effective interactions between two 
polymers of various lengths are plotted in Fig. g; the dis- 
tance r is scaled with the measured radius of gyration 
Rg. As expected, V2{r) has a Gaussian shape centered 
on r = 0. The potentials are almost indistinguishable for 
r/Rg > 1, but for smaller r the potentials differ slightly 
for different L. This is most pronounced at full overlap 
of the polymers, where V2{r — 0) decreases with length 
L. In the scaling limit L ^ oo, V2{r = 0) is expected to 
reach a finite value while for finite L, we expect V2{r = 0) 
to scale as: 

V2{r = 0) ex Lln(l - ac^/'^^"'^^) - Lln(l - -) (5) 

where a is a (negative) constant and the logarithmic 
term arises because P{r) scales linearly. This finite- 
size scaling behavior is confirmed in Fig.|3| and in the 
L — > oo limit this equation goes over to Eq. (^. Using 
a non-linear fit of the MC data to Eq. (J^) we estimate 
V2{r = 0) = 1.80 ± 0.05, a value slightly higher than the 
best ©(e^) RG calculations which give V2{r = 0) = 1.53. 
The difference is most likely due to a lack of convergence 
of the e expansion p^ . 
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FIG. 3. Finite-size scaling for the interaction between 
two isolated polymers at full overlap: V2{r = 0). By plotting 
exp [v2{r — 0)/L]v.s. 1/L, the agreement between the scaling 
relation of Eq. (H) and the simulations is demonstrated. 

The quantity r'^V2{r) is actually more relevant for 
the thermodynamic properties of polymer solutions than 
V2{r = 0) Q , and, as demonstrated in the inset of Fig. ||, 
the former varies less with L than the latter, such that 
for r'^V2{r) the scaling limit appears to be reached even 
for chains as short as L = 500. 

Similarly, with the effective pair-potentials we can cal- 
culate the 2nd osmotic virial coefficient: 



B2 = -2ti 



B2/ R^g should be independent of L in the scaling limit. 
As demonstrated in Fig.^, the scaling limit appears to 
be practically reached for L — 500. We estimate that for 
L -^ 00, B2/R^„ ~ 5.85 ± 0.05, which is consistent with 
other results obtained from simulations (B2/-R0 ~ 5.50 
||l|) or RG calculations (B2/i?| ~ 5.99 |42)). Note that 
although B2 scales as B2 ~ i??, as required by scaling 
theory H0,0|, this does not imply that the polymer- 
polymer interaction is hard-sphere like, as is sometimes 
implied in the literature p3|. 
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FIG. 4. The reduced osmotic virial coefficient B2/Rg v.s. 
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IV. EFFECTIVE POTENTIALS: POLYMER 
SOLUTIONS 



A. Deriving effective potentials from g{r) 

Having derived the effective potential between two iso- 
lated polymers, we now attempt the same for polymers 
in solution at finite concentration. Whereas for simple 
fluids, the interaction potential is generally independent 
of the thermodynamic state, this is not true for effective 
potentials in complex fluids. The latter typically follow 
from a coarse-graining procedure, which amounts to av- 
eraging out certain degrees of freedom, the individual 
microscopic polymer segments in the present case. The 
effective total interaction potential energy V^r({ri}; p^) is 
in fact a free energy which depends here on the polymer 
density pb — N/V and on the configuration {vi} of the 
polymer CM's. The bare pair-interaction term V2{r) can 
be defined as the effective potential between two isolated 
polymers, the bare triplet interaction term W3(ri,rj,rfc) 
can be defined for three isolated polymers and so forth. 
One could in principle calculate higher and higher or- 
der n-body terms, but this rapidly becomes intractable. 



Even if explicit expressions for each of the terms were 
obtained, the total interaction energy would be very dif- 
ficult to evaluate because the number of n-tuple coordi- 
nates increases exponentially. 

Instead, we follow a different route and approximate 
the pair and higher order terms by an effective, (state 
dependent) pair interaction v(r; pb) which is constructed 
to exactly reproduce the two-body correlations of the full 
underlying many-body system. In fact, it can be proven 
that for any given pair distribution function g(r) and 
density pb, there exists a corresponding unique two-body 
pair potential v(r; pb) which reproduces g{r) irrespective 
of the underlying many-body interactions in the system 
Q. Of course, g{r) will contain contributions not only 
from the bare pair-potential V2{r), but also from the three 
and more body terms. As a consequence, the effective 
pair interaction v{r; pb) will also be state dependent (in 
the polymer case, density dependent) and a new effective 
potential must be calculated for each density. Neverthe- 
less, the effective potential leads back to the true ther- 
modynamics of the full many-body system through the 
compressibility relation: 
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where h{k) is the Fourier transform (FT) of the pair cor- 
relation function h{r) — g{r) — 1, and c{k) is the FT 
of the direct correlation function. Using a variational 
argument, Reatto |^ has shown that v{r; pb) may also 
be viewed as the "best" pair representation of the true 
interactions. However this inversion approach says noth- 
ing about a possible volume term Vo{pb), in the coarse- 
grained total potential energy, which contributes to the 
e.o.s. , but not directly to the pair-correlations [Q. Of 
course the volume terms may still contribute indirectly, 
for example when they induce phase-transitions. 

The inversion of g{r) to extract v{r;pb) is a well 
known procedure and has been studied extensively in 
the field of simple fluids ^,0. We invert g{r) using 
the hypernetted-chain (HNC) closure, 



g{r) — exp{—v{r) + g{r) — c{r) — 1), 



(8) 



of the Ornstein-Zernike equation WSl . While the simple 
HNC inversion procedure would be inadequate for dense 
fluids of hard core particles, where more sophisticated 
closures or iterative procedures are required Bq,B^, we 
are able to demonstrate the consistency of the HNC in- 
version in the present case. 

We performed Monte Carlo simulations of A^ SAW 
polymers of length L = 500 on a cubic lattice of size 
M = 240 X 240 X 240. The number of polymers was varied 
from A^ == 100 {pb/p* = 0.54) to A = 6400 (pb/p* = 8.7). 
Note that at the highest density the monomer packing 
fraction is c « 0.23, meaning that the conditions for the 



semi-dilute regime, namely pb > p* and c « begin 
to be violated. At even higher densities the system will 
approach the melt regime where monomer packing ef- 
fects become important [^. More generally, for finite 
length SAW polymers, there is a limited density regime 
for which both conditions for the semi-dilute regime can 
be simultaneously satisfied. We find empirically that 
Rg ~ 0.39L°'^ for SAW polymers on a simple cubic lat- 
tice, so that the monomer packing fraction at the overlap 
concentration is given by 



4/L° 



(9) 



Thus, for L = 100 chains we find c* « 0.1 so that there 
is only a very small density range which might be called 
semi-dilute, while for L = 500 chains c* w 0.027 and 
a meaningful semi-dilute regime exists. The literature 
contains several claims of semi-dilute scaling behavior for 
SAW lattice polymers with L < 100, but, as the analysis 
above shows, these polymers do not have a semi-dilute 
regime large enough to derive scaling relations. An ex- 
ample of this is shown in Fig. H. 
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FIG. 5. This figure compares the e.o.s. Z — Ilb/pb from 
simulations for L = 40 SAW polymers |50|, wit h the e.o.s. 
for L = 500 SAW polymers (see sectiori |[ V B|) , and with 
the Bawendi- Freed (BF) e.o.s. for lattice models [^ll, which 
is accurate for larger values of c. The latter gives an in- 
dication where finite c correlations become important, and 
where one would expect the melt-regime to start for L = 40, 
L — 100, and L = 500 SAW polymers. In the melt regime 
the des Cloizeaux scaling law Z ~ {pb/p*y will break 

down M. Clearly, the L = 40 data does not follow the des 
Cloizeaux scaling law, demonstrating that there is no mean- 
ingful semi-dilute regime for L = 40 polymers, whereas there 
is one for L — 500 polymers. 



In the course of the simulations the CM of each poly- 
mer was tracked in order to construct the CM radial pair 
distribution function g{r). The latter is only known up 
to a cutoff radius r^, which corresponds to half the size 
of the simulation box (lattice size). For the inversion. 



we need g{r) for all r, so we employ the following iter- 
ative scheme to extend g(r). First we set g{r) = 1 for 
r > re and calculate the corresponding v{r; pb) by inver- 
sion. We then set u(r; pb) = for r > Vc and determine 
the corresponding g{r) for < r < oo by a regular HNC 
calculation, using a simple iterative procedure. The g{r) 
for r < re is then replaced by the measured g{r), and the 
process is repeated until convergence. For low density, 
g{r) and v{r; p^) converge very quickly, but for higher 
densities, say Pb/p* > 1, the convergence is slower, and 
the mixing factor of the old solution into the new one has 
to be increased to a value as large as 99%. In fact, be- 
cause of the finite box-size, the inversion process is under- 
determined, and our ansatz that v{r; pb) = for r > re is 
needed to find a unique solution. This is not unreason- 
able since we don't expect the interactions between the 
polymer coils to be significant beyond a distance a few 
times the radius of gyration. However, to make sure that 
this is actually the case, we found that relatively large 
simulation boxes were needed, with a lattice size of up to 
10 — 15i?g. This is especially important at high density, 
where the inverted potential becomes longer ranged and 
more sensitive to small changes in the radial distribution 
function g{r). In all our inversions we checked explic- 
itly that v{r\pb) becomes effectively zero for an r < re, 
confirming our initial ansatz. 

The resulting radial distribution functions g{r) are 
shown in Fig. 0, and are similar in shape to those of 
the pure Gaussian core model |^. As the density goes 
up the correlation hole at small r decreases in range and 
height. Except for a small maximum around r « 2i?g, 
the pair correlation functions do not show oscillations 
within the statistical noise of about 0.1%, for any den- 
sity considered here. 
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FIG. 6. The polymer CM pair distribution function g{r) 
calculated for L=500 SAW polymers and used to generate 
v{r;pi,). The x-axis denotes r/Rg, where Rg is the radius of 
gyration of an isolated SAW polymer. 




FIG. 7. The effective polymer CM pair potential v{r; pt) 
derived from an HNC inversion of g{r) for different densities. 
The X-axis denotes r/Rg , where -Rg is the radius of gyration 
of an isolated SAW polymer. Inset: The value of the effective 
polymer CM pair potential at r = 0, as a function of den- 
sity pb/p* ■ The maximum of the potential initially increases 
before decreasing at high concentration. 



The effective polymer-polymer potentials v(r;pb), ob- 
tained from the g(r)'s, are shown in Fig.M. Careful in- 
spection of the figure reveals that the effective pair po- 
tential is not very sensitive to the polymer concentration. 
The value at r = first increases slightly with pb, before 
decreasing again at the highest concentrations, as is de- 
picted in the inset of Fig H, while the range of v{r\ pb) 
increases with pb. A more subtle feature, highlighted in 
Fig. H, is that the effective potential becomes slightly neg- 
ative {O{\0-^kBT)) for r/Rg > 3 at the higher concen- 
trations. These effects become apparent only when large 
enough box-sizes are used. Although the negative tails 
seem very small, they are nevertheless significant since 
the thermodynamics depend on the integral of r'^v{r; pb). 
For example, leaving them out can easily induce a 5% 
change in the pressure. It is, therefore, paramount to 
include these effects in (quasi)-analytical representations 
of the effective potentials. For that reason, a simple fit to 
a Gaussian or a sum of Gaussians is not accurate enough 
to reproduce the structure and the thermodynamics of 
the SAW polymer systems and consequently, we chose to 
use an interpolation spline fit to describe the potentials. 
First, the raw effective potential data were fitted to a 
Gaussian: 



Vest{r) = floe 



(10) 



Subsequently, the difference /S.v{r; pb) = v(r; pb) — Vcst{f) 
was fitted by employing a least squares spline proce- 
dure with 8 nodes (the "dfc" routine of the slatec li- 
brary [p2|). The values of the nodes are not known in 



advance, except for the boundaries r = and r = re- 
Additional constraints on the sphne fit were: v{rc) = 0, 
dv{r = rc)/dr — and v{r — 0)/dr — 0. We optimized 
the sphne fit by moving the nodes on the x-axis using a 
Monte Carlo procedure. The parameters for the fits are 
available elsewhere lls^. 
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FIG. 8. The negative part of the effective polymer CM 
pair potential v{r; pt.) derived from an HNC inversion of g{r) 
for different densities. The x-axis denotes r/Rg, where Rg is 
the radius of gyration of an isolated SAW polymer. 



Note that in Fig. M, the polymer-polymer interaction 
v{r;pf,) is plotted v.s. r/Rg, where Rg is the radius 
of gyration of isolated polymers in the infinitely dilute 
limit. In a dense solution, the effective radius of gyration 
of the polymers contracts according to the power-law, 



Rn 



Pb 



1/8 



|||,||l, as shown in Fig. |. 




FIG. 9. The effective radius of gyration for L = 500 SAW 
polymers decreases as a function of density pb/ p* ■ Inset: At 
high densities the effective radius of gyration asymptotically 



follows the scaling law Rg 



^1/8 



The accuracy of the effective potentials are tested by 
performing a direct Molecular Dynamics (MD) simula- 
tion of the "soft colloids" interacting via v{r;pb). In 
Fig. no the pair distribution function gMBif) from MD 
simulations is compared with the original SAW g(r) for 
two densities in the semi-dilute regime. The difference 
between the two distribution functions shows an oscilla- 
tion at small r. Because this occurs in the same way for 
both densities it is possibly introduced by the inversion 
procedure. Even so, the difference between the two dis- 
tribution functions is still typically less then ±0.01. We 
conclude that the HNC inversion procedure yields very 
accurate effective potentials for soft particles, capable of 
describing the structure of the fluid with an absolute er- 
ror of less than ±0.01. 




FIG. 10. The g{r) of a system interacting via the effective 
potential v{r; pb) compared with the CM pair distribution of 
a SAW simulation for two polymer concentrations. The dif- 
ferences are shown in the lower panel and are typically less 
than ±0.01. 



In a previous paper we have shown that the HNC clo- 
sure is very accurate when applied to the Gaussian model 
p5[ , whereby particles interact via the repulsive poten- 
tial v{r) = eexp[— a (r/i?g) ], and is in fact quasi-exact 
in the regime relevant to the effective potentials shown 
in Fig. R Eolpfl]. Even the much cruder RPA closure, 
c{r) = —v(r), yields semi-quantitatively accurate results 
for correlations and thermodynamics in the regime of in- 
terest. Thus polymer solutions in the dilute or semi- 
dilute regime fall into the class of mean field fiuids ac- 
cording to the nomenclature introduced in Ref. flOl] . 

The inversion procedure guarantees that the two-body 
correlations are accurately reproduced by the effective 
potential, but this does not necessarily imply that higher 
order correlations are also well represented. As a first test 
we performed preliminary simulations of the three-body 
bond-order correlation functions gj,{r,9,Lp) for both full 
SAW walks and our soft particles. The two approaches 



lead to identical results within statistical errors, imply- 
ing that higher order correlations are much more accu- 
rately reproduced than one might initially expect. We 
have also performed some preliminary calculations of the 
three-body interaction W3(ri, r2, ra). Even at full overlap 
of the three centers of mass, the three-body interaction 
term is only about 10% of the pairwise interaction. This 
is consistent with the results found for star-polymers [^ , 
and was fore-shadowed by the relatively weak density de- 
pendence of the effective pair interaction v(r;pb). 

Besides accurately describing the structure, it is also 
important that the thermodynamics are captured by the 
effective potential. In the next section we therefore focus 
on the equation of state (e.o.s.) for polymer solutions. 



B. Equation of state 

1. Equation of state from direct SAW simulations 

We measured the e.o.s., Ilb/pb, directly for a SAW 
simulation by using the thermodynamic integration ap- 
proach of Dickman pM]. In this method the bulk (os- 
motic) pressure lib is measured by taking the derivative 
of the free energy F with respect to volume of a system 
of SAW polymers between two hard walls. The polymers 
live on a rectangular cubic lattice of size M = HxDxD, 
which is periodic in the y and z directions. The two 
walls are represented by an infinitely repulsive poten- 
tial at a; = and at x = H + I, so that the polymer 
segments cannot penetrate the walls. The volume of a 
lattice can only change discretely, and the free energy 
derivative changes to a finite difference 

d\nZ{N,L,D,H) ^^d\nZ{N,L,D,H) 

lib = TT-rz = J^ 



DM 



dH 



w D-^ (In Z{N, L, D, H) - In Z{N, L,D,H- 1)) . (11) 

The model is modified by associating an additional repul- 
sive potential — In A with each occupied site in the plane 
X — H, where < A < 1. The partition function then 
becomes 



Z{N,L,D,HX)= Y^ e 

polymer conf 



-U \n 



•A"«, 



(12) 



where nn = D^pni^) is the number of occupied sites in 
the X = H plane, and PhW is the corresponding number 
density in this plane. The pressure can now be estimated 
as 



Tl, = D 



-^'."-(^ 



\x^ (13) 
A 



We performed SAW simulations of polymers with length 
L = 500 on a M = 160 x 100 x 100 cubic lattice for 
iV = 50, 100, 200, 400, 600 and 800. For each density we 



determined the value of PhW for 5 different values of 
A, corresponding to the abscissae of a 5 point Gaussian 
quadrature which was used to evaluate the integral in 
Eq. ( |l3| ) The resulting e.o.s. is plotted in Fig. |ll|. 
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FIG. 11. Log-log plot of the e.o.s. Z — U/pb as a function 
of the density for L — 500 polymers. The soft-particle e.o.s. 
gives a good representation of the full SAW polymer simula- 
tions. At the highest densities there is a slight deviation from 
the expected des Cloizeaux {pb/ p*Y scaling law which 

we attribute to the effects of a finite monomer concentration 
c. Also shown is the RG e.o.s. of Ohta and Oono |6d,pl|. 



2. Equation of state from the soft-particle picture 

To calculate the e.o.s. within the soft-particle picture, 
we use the compressibility relation (|7|), which must now 
be integrated w.r.t. the density: 



fPb 

Tib{pb)= / {l-p'c{Q,p'))dp'. 
Jo 



(14) 



We used the quasi-exact HNC approximation to calcu- 
late c(r) from the inverted effective potential v{r; pb) for 
several state-points, fitted the values of c(0;pb), and in- 
tegrated w.r.t. density. As demonstrated in Fig. |ri|, the 
e.o.s. is very close to the one obtained by direct SAW 
simulations, immediately suggesting that our inversion 
procedure indeed reproduces the true thermodynamics 
of the full many-body system. This success also implies 
that the volume terms are small, possibly smaller than 
the statistical error in our present simulations and in- 
verted potentials. In fact, using a simple scaling theory, 
Likos has argued that the contribution of volume terms 
to the e.o.s. scales as {pb/ P*)^^^ in the semi-dilute regime, 
and so contributes little to the full e.o.s. [B9[. 

Also shown in Fig. ni^ is the RG result by Ohta and 
Oono pO| ; we use a slight improvement with correct ex- 
ponents [ pl[ . The one remaining fit parameter is de- 
termined by the second osmotic virial coefficient B2 for 



L = 500 SAW polymers, a procedure similar to that used 
when comparing to experiment p% . The agreement is 
seen to be fairly good, although the SAW e.o.s. is some- 
what higher than the RG results. This is most likely due 
to the fact that the monomer density c is not zero, which 
induces small corrections to the full scaling limit (see the 
discussion in Section IV A). 
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FIG. 12. Linear plot of the the e.o.s. Z = Ht/pb as a func- 
tion of the density for L = 500 polymers. Several approxima- 
tions to the e.o.s., discussed in the text, are compared. 



density derivative terms in the virial equation, we over- 
estimate the e.o.s. compared to the full SAW simulation. 

The density dependence can be neglected even further 
by simply taking the pb ^ form of the pair potential, 
V2{r), and applying it at all densities. The resulting e.o.s. 
now underestimates the e.o.s. when compared to the full 
SAW simulation, as demonstrated in Fig. |l2|. We note 
that a similar approach was employed in recent work on 
the phase-behavior of star-polymers, where the pb ^ 
limit of the pair potential was used to calculate the struc- 
ture and phase-behavior at finite concentration [Q . 

Finally, we comment on the common practice of ex- 
tracting osmotic virial coefficients from the measured ex- 
perimental e.o.s. Firstly, the virial equation has a very 
small radius of convergence for soft-core fluids |4^] . Sec- 
ondly, the range of the effective potential v{r; pb) in- 
creases with density. These two effects imply that a naive 
linear fit to all but the very lowest polymer densities will 
lead to an overestimate of the true osmotic second virial 
coefficient i?2. 



V. EFFECTIVE WALL-POLYMER POTENTIALS 



A. Polymer coils near a wall 



Instead of the compressibility route, one could also use 
the virial route to the e.o.s. ifisi: 



Hfc ^^ ^ dVojpb) 
Pb dpb 



1 



N 






dv{rij]pb) 
dri. 



-3pj 



dvjrij^pb) 
dpb 



(15) 



which includes not only the density dependence of the ef- 
fective pair potential, but also the density dependence of 
the volume term. The full density dependence of the po- 
tentials is at present hard to calculate, so instead we ini- 
tially ignore the density derivative and the volume terms. 
First, we directly measured the e.o.s. of the soft-particle 
fluid by a MD simulation with the spline-flt potentials. 
The pressure follows from the usual virial theorem, when 
the density derivatives in Eq. ( [l5|) are neglected. The 
e.o.s. from this approach is depicted in Fig. n2, and com- 
pared to the simple mean-field (MF) form: 



Zmf = 1 + -v{0;pb)p, 



(16) 



which gives a good fit to the simulations, as expected for 
soft-core fluids in the MFF regime Eo]. Here v{0; pb) is 
the fc = component of the FT of the pair interaction. 
However, by including only the explicit density depen- 
dence of the effective pair potentials while ignoring the 



Polymer coils near a non-adsorbing hard wall exhibit 
a depletion layer due to entropic effects. This is true 
even for ideal Gaussian polymers, and if one were to 
model these by effective CM potentials, the polymer- 
polymer potential would be zero, but there would still be 
a polymer- wall potential of the form (j){z) = ln(p(z)/pf,), 
where p{z) is the CM density profile near the wall and 
Pb is the uniform density far from the wall (see the Ap- 
pendix for more details). Thus a complete description of 
polymer coils in confined geometries requires not only the 
polymer-polymer interactions derived in the previous sec- 
tion, but also effective polymer- wall potentials (f){z;pb). 

We follow a strategy similar to that used in the homo- 
geneous case, and first calculate the wall-polymer den- 
sity profile p{z), from which we then extract an effective 
potential (f){z;pb). Using the same explicit SAW poly- 
mer model as in Section ||, we performed MC simu- 
lations of polymers of length L = 500 on a lattice of 
size M = 160 X 100 X 100 with hard walls at x = 
and X = 160. The polymer segments were not allowed 
to penetrate the walls. The simulations were done for 
N — 50, 100, 200 and 500. During each simulation, we 
computed the density profiles p{z), where z denotes the 
distance of the polymer CM from the wall and pb is the 
bulk density far from the wall. The normalized profiles 
h{z) = p{z)/pb—l, for different bulk concentrations pb/p* 
are shown in Fig. O. 



/ V 



N 



-0.5 




Pb/P* = 
p^/p* = 0.67 
p,/p* = 1.28 
p^/p* = 2.49 
p^/p* = 6.05 



12 3 

FIG. 13. The wall-polymer CM density profile 

h{z) = p{z)/pb — 1 for SAW polymers at different bulk con- 
centrations. From h(z) we can calculate the corresponding 
polymer adsorptions F and find — F — 0, 0.094, 0.13, 0.16, and 
0.20 in units of R~^ respectively. The relative adsorptions 
are -T/pt, = 0.84, 0.59, 0.41, 0.27, and , 0.14 respectively, and 
decrease with increasing density as expected. 



The polymer coil adsorption T is defined by: 



d^ 



= Pb 



h{z)dz, 



(17) 



where ft'^^ /A is the excess grand potential per unit area 
and /i the chemical potential of the polymers. As the 
density increases, more polymer is adsorbed at the wall 
as expected, but the relative adsorption, F/pb, decreases. 
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FIG. 14. 

The wall-monomer density profile </)(z; pt) = p{z)/pb — 1 for 
the same set of densities as in Fig.y_3| Inset: A magnifica- 
tion of the region where there is a small correlational bump 
in the density profiles. The height is less than 1% of the total 
density while the range is about Rg , implying that the bump 
arises from polymer-polymer correlations. 



The normalized monomer density profiles for SAW's 
are shown in Fig. n4 for the same polymer densities as 
the CM profiles shown in Fig. |3[ As expected, the profile 
moves closer to the wall for higher density; the width of 
the monomer depletion layer shifts from around Rg at the 
lowest densities, down to values dictated by the segment 
correlation length |2] in the semi-dilute regime. Although 
the profiles do not show such a clear correlation-induced 
oscillation as the CM profiles, there is nevertheless still 
a small maximum in the depletion layer as illustrated in 
the inset of Fig n4. The peak in the monomer profile is 
less than 1% of the bulk density, and seems to decrease 
for higher overall polymer concentration. The range is 
about Rg, implying that it arises from correlations be- 
tween polymer coils. We observe only one peak, although 
due to statistical noise, we cannot rule out the possibility 
of more oscillations in the density profiles. Recently self- 
consistent field calculations, valid for polymers in a theta 
solvent, found a similar small oscillation in the monomer 
profiles il. 



B. Deriving <j)(z; pi,) from p(z) 



From a knowledge of the concentration profile p{z), 
and the bulk direct correlation function between polymer 
CM's, ci,{r), one may extract an effective wall-polymer 
potential (f){z; pb) by combining the wall-polymer OZ re- 
lations iQ with the HNC closure. For a binary mixture 
of two components labeled and 1 , in which component 
is infinitely dilute (xq — > 0), the Ornstein-Zernike equa- 
tions become lEsI 



/iii(l, 2) = cii(l, 2) + pb /iii(l, 3)cn(2, 3)d3 (18a) 

hw{l, 2) = cio(l, 2) + pb I /iii(l, 3)cio(2, 3)d3 (18b) 

/ioi(l,2)==coi(l,2) + pb /"/ioi(l,3)cn(2,3)d3 (18c) 

/ioo(l,2) = coo(l,2) + p6 //ioi(l,3)cio(2,3)d3 (18d) 



In the limit Rq —f cx3, Eq. (18c) becomes an equation for 
the wall-density profile, sometimes called the wall-OZ re- 
lation: 

h{z) = coi(z) + pb f dr'hoi{z')ci,{\r - r'|), (19) 



where h(z) = p{z)/pi, — 1. The wall-OZ equation can 
be solved, given the bulk correlation function Cb{r), and 
a closure relation. In Section. IV A we showed that the 



HNC closure gives excellent results for effective polymer- 
polymer interactions, and it is therefore natural to apply 
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the same approximation here. Combining Eq. (@) with 
Eq. ( pj| ) we obtain 

<P{z-pb)^^^'^{z-pb)+pb fdr'h{z')cb{\r-r'\). (20) 
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FIG. 15. The wall-polymer potential (j)(z; pb) as obtained 
Irom the inversion ol h(z) via the HNC expression, Eq. (EOl). 



The first term is the usual potential of mean force 
(/>*^^(z; Pf,) = ^ In [/o(2:)//Ob], to which (j){z;pi,) would re- 
duce in the pf, ^ limit, while the second term arises 
from correlations between the polymer coils next to the 
wall. An identical equation results from the HNC den- 
sity functional theory (DFT) approach [ p5[ , and a simi- 
lar one, with Cb{r) replaced by v{r; pb) obtains if a mean 
field DFT is used. In contrast to simple fluids, where 
Eq. (^0|) is not very reliable, the wall-HNC closure works 
remarkably well for the Gaussian core-fluid in the regime 
relevant to polymer solutions |Q. Using the Cf,(r) ex- 
tracted from the earlier bulk simulations of ^(r) (see Sec- 
tion IV A ), together with Eq. (|20|), we are able to extract 
4>{z\ ph) from the density profiles. In order to calculate 
the integral in Eq. d2fl), we use the procedure outlined by 
Sullivan and Stell pq] . In contrast to the inversion of the 
bulk g{r), where we had to iterate until convergence, the 
wall-polymer inversion requires only one step since Cb{r) 
is given once and for all. Results for various bulk concen- 
trations are plotted in Fig. UB. The range of the effective 
wall-polymer repulsion increases with increasing concen- 
tration, while the density profiles actually move in closer 
to the wall. The compression and enhanced correlation 
in the density profiles with increasing density resembles 
that of the pure Gaussian core fluid in a fixed external 
potential Eo], but the effect is less pronounced in the 
former case since for polymer solutions the wall-polymer 
potential becomes more repulsive with density. This is 



due mainly to the correlation term, which is nearly lin- 
ear in Pb, and so becomes relatively more important as 
the density increases. Nevertheless at shorter distances 
the (j)^^ {z]pb) term still dominates. The importance 
of including both the potential of mean force, and the 
correlation-induced component of the effective potential 



is demonstrated in Fig, 16, At very low densities the po- 



tential of mean force is adequate, but at higher densities 
the correlation term becomes increasingly important. 
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FIG. 16. Comparison between the contributions to 

the effective wall-polymer potentials from the potential of 
mean force (solid lines) and from the correlation part 
(dashed lines) (cf. Eq. (EOl)) for polymer concentrations 
Pf,/p* = 0.67, 1.28, 2.49 and 6.05. From top to bottom the 
solid lines correspond to increasing density and the dashed 
lines correspond to decreasing density. The solid line with 
the small squares denotes the potential of mean force for in- 
finitely diluted systems. 



The effective potentials decay exponentially, and to ob- 
tain a useful analytic form for the effective potential, the 
logarithm of ^(z; pb) can be fitted to a cubic polynomial, 
which describes the potential very well. However, as in 
the bulk case, the wall-polymer potential (j){z] pb) has a 
small negative component that cannot be described by 
an exponential function. Although in this case the tail is 
probably not very important, in order to be consistent, 
we fit (/)(z; Pb) by a least squares spline fit similar to the 
one described Section. 



IV A. The parameters for this fit 



are available elsewhere fB3f 



C. Consistency of the wall-polymer inversion 

To test the validity of the inversion procedure for the 
wall-polymer p{z), we performed Molecular Dynamics 
simulations of a system of "soft colloidal" particles in- 
teracting with each other via the effective pair potential 
v{r, Pb) and with a wall via the inverted potential 4>{z; pb) 
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for the appropriate bulk concentration pb- Such effective 
potential simulations are at least an order of magnitude 
faster than simulations of the original SAW model. The 
resulting concentration profile of the effective particles is 
shown in Fig. 0^ for one density; it agrees to within an 
absolute error of roughly ±0.02 with the p{z)/ ph obtained 
from the detailed SAW simulations. The corresponding 
adsorption F also differs by less than 1% from the value 
obtained by the SAW simulation, thus demonstrating the 
adequacy of the soft colloid representation of the interact- 
ing polymer coils, and the accuracy of the HNC inversion 
for polymer coils near a hard wall. 




FIG. 17. The profile h{z) of a system of soft colloids near 
a wall (dashed line). The particles interact with each other 
via v{r), and with the wall via (jti^'^Pb)- This is compared 
with the wall-CM distribution of an explicit SAW simulation 
(solid line) . The difference is shown in the lower panel and is 
less than ±0.02. 



VI. DEPLETION POTENTIAL BETWEEN TWO 
WALLS 

A. Full SAW simulations 

One of the aims of this work is to show that the soft 
particle description of the polymers provides a useful 
route to the colloid-colloid depletion potential in mix- 
tures of colloidal particles and non-adsorbing polymers. 
Calculating these depletion interactions poses a severe 
test of the soft colloid representation. 

As a first step we calculate the depletion potential be- 
tween two planar walls, which can ultimately be applied 
to spherical colloids through the Derjaguin approxima- 
tion. We confined the polymers within a slit of width 
d, and, using direct grand-canonical simulations of the 
full SAW polymer model, we computed the osmotic pres- 
sure exerted by the polymer coils on the walls. The in- 
sertion of polymers was achieved by the configurational 



bias Monte Carlo technique ||3^]. The (osmotic) pressure 
n((i) was calculated for different values of the spacing d 
between the walls by a thermodynamic integration tech- 
nique similar to the one explained in Section fVB. De- 
tails of these simulations can be found in Ref. |)7[. The 
interaction free energy per unit area A, AF/A, is then 
obtained by integrating the osmotic pressure as a func- 
tion of d: 



AF{d)/A 



dz(n(z)-n(oo)). 



(21) 



where n(oo) denotes the bulk osmotic pressure lib- These 
explicit SAW simulations are rather computer intensive, 
and were only carried out for L = 100 [p7[. 



B. Effective potential simulations 

In the soft colloid picture, the interactions of the poly- 
mer CM's with each other, v{r;pb), and with a wall, 
(f>{z; Pb), are calculated once with the HNC inversion pro- 
cedures from the g{r) and p{z) of a full SAW polymer 
simulation at the bulk density pb- These potentials are 
then used in grand-canonical MC simulations of soft par- 
ticles between two walls. The imposed chemical poten- 
tial is chosen such that for infinitely separated walls the 
bulk density is recovered. The (osmotic) virial pressure 
is measured as a function of wall separation d, and the 
interaction free energy per unit area AF/A, is again ob- 
tained by integration of the pressure via Eq. ( ^ . 

In Fig. O the soft colloid depletion interaction is com- 
pared to that of the "exact" grand-canonical MC simula- 
tions oi L ~ 100 SAW polymers, for three different den- 
sities, pb/p* — 0.28, 0.58 and 0.95. The two approaches 
are in good agreement, but the soft colloid calculations 
are at least two orders of magnitude faster than the SAW 
simulations. As expected, the depth of the potential in- 
creases, whereas the range of the interaction decreases as 
the density increases |Q . At the two lowest densities, the 
two approaches agree very well, but for Pb/ p* = 0.95 they 
differ slightly around z = 2Rg where the soft particle pic- 
ture shows a larger repulsive barrier. The barrier height 
is, however, small compared to the attractive minimum 
at contact, which agrees well with the "exact" data, as 
does the slope of the attraction. 

Liquid state theories for fluids with repulsive particle- 
particle interactions predict a repulsive barrier p9| , so it 
is not surprising that the soft particle picture shows a 
small repulsive barrier as well. Instead, it is the lack of 
a significant barrier for the pure SAW polymer simula- 
tions which requires explanation. We trace this effect to 
the breakdown of the "potential overlap approximation" 
(POA) described in the Appendix. Under close confine- 
ment, the interaction of the soft particles with two par- 
allel walls a distance d apart can no longer be written as 
the sum of the two individual wall-particle interactions 
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as would be the case for simple liquids. This is caused 
mainly by the deformation of the polymers due to the two 
walls, and also holds for ideal polymers. The failure of 
the POA can be clearly seen in Fig. |l9|, where we compare 
the pressure (or force) profiles for the SAW calculations 
and the effective potentials. In the soft particle picture 
the pressure starts to rise at a larger inter-wall distance 
than the pressure for the SAW polymers, an effect also 
seen when non-interaction polymers are represented by 
an effective particle representation based on the CM (see 
the Appendix). Note that the over and under- estimates 
of the pressure cancel each other, so that the free energy 
at contact, AF(0), for the effective potentials is in good 
agreement with the SAW calculations. 
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FIG. 18. Depletion free-energy AF{d)/R^ between two 
plates separated by d, for three densities, pt/p* ~ 0.95, 
Pb/p* ~ 0.58, phi p* = 0.24. The symbols denote the "exact" 
MC simulations, while the dashed, dash-dotted and solid lines 
are the soft-colloid simulations for the same densities. Inset: 
Af (ii)/A_F(0) for the SAW simulations, the solid lines are 
to guide the eye. The long-dashed line is the ideal Gaussian 
polymer result calculated in the Appendix. Note that the 
range decreases with density, and that, even for the lowest 
density, the AO ideal polymer approximation overestimates 
the interaction range. 



1.25 



0.75 



2, 

CL 




^--©SAW MC simulations 

► Soft particle MC simulations 



0.25 



-0.25 



2 4 6 8 

d/R, 

FIG. 19. Comparison of normalized the depletion force 
per unit area (pressure) Pid^jW^ for SAW polymers and soft 
particles between two plates separated by d for a density 
Pb/p* = 0.95. 



The MC simulations for the soft colloid model were 
carried out with effective wall-polymer and polymer- 
polymer potentials appropriate for L = 100, since longer 
polymers are not easily handled in the full SAW model. 
However, we checked that the data obtained with effec- 
tive interactions appropriate for longer polymers (L — 
500), are very close to the L = 100 results, as is shown 
in Fig. |20. Therefore, we are confident that we are close 
enough to the scaling regime for the properties under 
consideration. 
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FIG. 20. Depletion free-energy /S.F{d)/E?g between two 
plates separated by d based on the soft particle representa- 
tion for polymers of length L = 100 and L — 500. Here 
Pb/p* = 0.95. 
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C. Comparison with the Asakura-Oosawa Approach 

The first (and still most popular) approach to the de- 
pletion interaction in coUoid-polynicr mixtures was pio- 
neered by Asakura and Oosawa in 1954 |25|, when they 
approximated the polymers as ideal (Gaussian) , and cal- 
culated the induced attraction between two walls. We 
shall refer to this neglect of polymer-polymer repulsion 
as the AO approximation, in contrast to the AO model, 
where a further step is taken and the polymers are ap- 
proximated as inter-penetrating spheres of radius Rg p4| . 

The exact depletion potential induced by ideal poly- 
mers between two plates of area A a distance d apart is 
given by: 



AF{d)/A^PbAV,d{d), 



(22) 



where AVid{d) is the gain in volume accessible to an ideal 
Gaussian polymer of size Rg , due to overlap of the exclu- 
sion volumes close to the plates. This can be exactly cal- 
culated as shown in the Appendix. To treat interacting 
polymers, a widely used phenomenological improvement 
(see for example ref . Q ) replaces the ideal polymer den- 
sity by the bulk osmotic pressure lib of the interacting 
polymers in the left over free-volume: 



AF{d)/A^UbAV^d{d). 



(23) 



In Fig. ^ we plot these two versions of the AO ap- 
proximation for the largest density considered above, 
Pb/ P* = 0.95, and compare them to the effective poten- 
tial and "exact" SAW simulation results. The two ap- 
proaches result in rather poor representations of both the 
depth and the range of the true potential, even though we 
are technically not yet into the semi-dilute regime where 
one might expect them to break down (see also the inset 
of Fig. |l8|). For the lower densities the AO approximation 
works somewhat better, as expected. 
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FIG. 21. Depletion free-energy AF{d)/R^ between two 
plates separated by d for pbj p* = 0.95 . Circles are the "ex- 
act" MC simulations of SAW polymers. The long-dashed and 
dash-dotted lines denote the two AO approximations men- 
tioned in the text. The short-dashed line denotes the more 
accurate wall-HNC approximation of Eq. (gj) , which is in fact 
very close to the the simulations in the soft particle picture 
shown in Fig. [isl 



D. HNC Wall- Wall Approximation 

Following arguments similar to those u sed t o derive 
the wall-polymer HNC equations of section VB , one can 
also derive an HNC type equation for the depletion in- 
teraction free energy per unit area between two walls 
separated by a distance d |7Q]: 



^^Pb I h{s)h{d — s)ds 

..Z — oo 

/•OO 

+ Pb / h{d - s) [0(s; Pb) ~ 0*^^ (s; pb)\ ds. (24) 



AF{d) 
A 



Here h{s) — p{s)/ pb — 1 is the single wall density pro- 
file, 0(s; Pb) is the corresponding effective wall-polymer 
potential, and (j)^^ {s; pb) is the corresponding potential 
of mean force. The first term on the r.h.s. of Eq. ( [24| ) 
represents the density overlap approximation discussed 
in the Appendix, and is the only contribution in the case 
of ideal Gaussian polymer coils. The second term arises 
from correlations between the polymer coils, and domi- 
nates the first term for larger densities. Note that only 
information from one single wall enters into this HNC 
wall-wall approach. We use the effective wall-polymer po- 
tential (j){z; Pb) and the related density profile h{z) from 
the soft-colloid picture together with Eq. (p3) to derive 
the HNC wall-wall depletion free energy. As shown in 
Fig Eh, this compares well with the MC simulations of the 
soft particles except at short distances where a small de- 
viation develops that can be traced to the fact that only 
information from a single wall is used. A more promising 
approach, without this shortcoming, would be to directly 
use the MF or HNC DFT approaches applied in ref. 0] 
to Gaussian-core potentials. 



VII. DISCUSSION AND CONCLUSION 

The coarse-grained representation of polymer coils as 
soft colloids, put forward in this paper, has proved 
very reliable. The effective polymer-polymer and wall- 
polymer interactions obtained by a systematic inversion 
procedure based on fluid integral equations, yield pair 
distribution functions and concentration profiles which 
agree closely with the results from simulation of the full 
SAW segment model, while allowing a massive reduction 
in computer time compared to the lattice simulations. 
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Much of the success of the present coarse-graining pro- 
cedure hes in our finding that the optimum effective pair 
potential between the CM's of neighboring coils does not 
depend strongly on polymer concentration, and is rea- 
sonably close to its infinite dilution limit. The effective 
polymer-polymer and wall-polymer interaction lead to a 
rather accurate description of the depletion interaction 
between two hard walls, despite the implicit potential 
superposition assumption and the fact that the coarse- 
graining procedure in its present form does not allow 
for the deformation of the polymer coils, away from the 
spherical shape, in the vicinity of an impenetrable sur- 
face. Such shape fluctuations are allowed in the alter- 
native procedure by Murat and Kremer |g^ but the re- 
markable agreement between the original full SAW model 
and the coarse-grained model illustrated in Fig. |lj seems 
to indicate that shape deformation of confined polymers 
may not be a crucial factor to reproduce concentration 
profiles. 

The present inversion procedure yields concentration 
dependent effective pair potentials, but does not pro- 
vide direct access to the internal free energy of polymer 
coils [g8[ , which plays a role rather similar to that of the 
"self energy" or volume term of electric double-layers in 
charge-stabilized colloidal dispersions or solutions of star 
polymers Eq]. This concentration-dependent term con- 
tributes to the osmotic equation of state, but the good 
agreement between full SAW model simulations and the 
results based on the effective pair interaction without the 
volume term would indicate that the concentration de- 
pendence of the latter is weak. 

Finally, it must be stressed that the present inversion 
procedure is by no means restricted to the simple SAW 
model of non-intersecting polymers. We are in fact plan- 
ning to extend the coarse-graining procedure to the case 
where the segment-segment coupling has an attractive 
component to describe the situation of polymer coils in 
poor solvent. The case of semi-dilute solutions of poly- 
mers of different lengths will also be considered within the 
same theoretical framework with the objective of study- 
ing possible demixing, as suggested by our recent inves- 
tigation of binary Gaussian-core systems l^^l- -^ final 
extension is to consider explicitly colloid-polymer mix- 
tures, by determining the effective hard-sphere/polymer 
potential along the lines set out in this paper. The gen- 
eral methodology should, more generally, be applicable 
to dilute and semi-dilute solutions of linear, branched or 
star polymers in confined geometries. 
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APPENDIX A: DEPLETION POTENTIAL FOR 
IDEAL POLYMERS 

In this Appendix we pursue a programme similar to 
that of the main text, but now for the simpler case of 
ideal Gaussian polymer coils of size Rg. Consider two 

apart. In 



LxLy a distance L^ 



parallel walls of area A 

the limit Lx,Ly » Rg, the full partition function for a 

single polymer is given by |25[| : 



Zi — L^ 



LyLz 



E 

0=1,3,. 



p^ 



■ cxp 



TT^Ry 



(Al) 



From this, various properties, such as the depiction in- 
teraction between two walls, can be exactly calculated. 

Similarly, from the underlying Green's function (see 
e.g. p. 19 of ref. [pj), the polymer end-point and mid- 
point density distributions near a single wall are found 
to be: 



Pind(^) 






Pb 



= erf(^ 

\2Rg 



= erf 



Pb 



\V2Rn 



(A2) 
(A3) 



As shown in Fig. E3, both show a clear depletion layer 
where polymer configurations are entropically excluded 
near the wall. In principle the CM distribution could also 
be calculated, but we have not yet succeeded in finding 
an analytic expression for it. Instead a polymer lattice 
model simulation with L = 500 was used to generate the 



Pcm(-^) depicted in Fig. 
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FIG. 22. Normalized densities p^^'{z)/pb for ideal poly- 
mers near a wall. Also shown is the normalized monomer 
density |7l| which falls, not unexpectedly, between the end 
and mid-point densities. The CM profile is the steepest since 
the CM cannot easily approach the wall as closely as the end 
or midpoints can. 



Using either the end-points, mid-points or the CM as 
the centers of "effective particles", an analogy with the 
full "soft particle" picture can be made. The particle- 
particle effective interaction is of course zero, while the 
effective particle-wall interaction can be inverted from 



the densities in Eqs. A2 and A3 with the result 



CW = -ln 



erf 



*^LW--21n 



a(i) 



rJ\ 




z 


] 


V2R 


J\ 



(A4) 
(A5) 



The potential (/)J^J^,^(z) for the CM profile can be obtained 
numerically. In contrast to the soft-particle picture for 
interacting polymers, the inversion here is trivial, since 
for ideal particles the effective potential is simply the po- 
tential of mean force. The interaction 0^|^(z) for the 
Gaussian particles is in fact quite similar to the poten- 
tial of mean force for SAW polymers at infinite dilution, 
shown in Fig. Qq. But, whereas the wall-polymer inter- 
action for interacting polymers changes with the bulk 
density, the (/)^^\z) for Gaussian chains is independent of 
density. 

Within the effective particle picture, two ways of cal- 
culating the interaction between two parallel walls are: 
(1) the Potential Overlap Approximation (POA): 
Here the partition function of the effective particles con- 
fined between two walls, a distance Lz apart, is calculated 
for a total external potential given by: 



(z)-0(i)(z) + (/.«(L,-z) 



(A6) 



where z is the distance from one of the walls. For simple 
atomic or molecular fluids this superposition approxima- 
tion would be exact and lead to the correct partition 
function and related equilibrium properties. However, 
for effective particles this is not necessarily the case as 
we shall see later on. 

(2) the Density Overlap Approximation (DOA): 
Here the density between two parallel walls is approx- 
imated by the product of the densities near a single wall: 



p,piz)=p'^'>{z)p^''{Lz-z) 



W( 



(A7) 



In contrast to the POA, this approximation is incorrect 
even for simple atomic or molecular systems, although it 
is sometimes a useful first approximation. On the other 
hand, for ideal particles the POA and DOA approxima- 
tions are equivalent. 



The original Asakura Oosawa model [Q approximates 
the density profile p'^^^ (z) next to a single wall by a step 
function of range R — Rg. The depletion potential is 
then calculated within the DOA. This can be improved 
by adjusting the range of the step-function such that it 
excludes exactly the same amount of polymer as the true 
density profile. For flat walls this implies a step-function 
of range R = 2Rg 
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FIG. 23. A comparison of the normalized depletion inter- 
action /S.F{Lz)/Api, as a function of the wall-wall separation 
Lz/Rg for various approximations discussed in the text. The 
X-axis denotes. 



In Fig. g^ we plot the depletion free-energy induced by 
non-interacting polymers between two walls a distance 
Lz apart. The exact expression was first calculated by 
Asakura and Oosawa pq] ; here we approximate it by the 
following simple analytical expression: 



^^(^.). Jj.-Ji-V'* ;i.<"32«. 



A 

AFjLz 
A 



= 0; 



Lz > 4.332i?„, 



(A8) 



which arises from taking only the "ground state" of the 
partition function in Eq. (Al), and cutting the potential 
off where it crosses 0. This approximation is so accu- 
rate that the difference with the exact interaction cannot 
be resolved in Fig. |2^. The effective particle representa- 
tions, based on end-points, mid-points or the CM, pro- 
vide a fairly good approximation to the full depletion 
interaction, while the two versions of the AO model do 
not perform as well. 
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exact pressure profile 
end-point effective particles 
mid-point effective particles 
CM effective particles 



L/R„ 



FIG. 24. Comparison of the reduced pressure P{Lz ; i-i)/Ii-b 
between the walls as a function of wall- wall separation L^/Rg. 
Hb is the bulk pressure. 



The pressure profiles shown in Fig. |2j demonstrate 
that the end and mid-point effective particle approaches 
underestimate the steepness of the pressure profile, while 
their CM counterpart overestimates the steepness. The 
differences between the effective particle representations 
and the exact results arise from our use of the POA. 
Because the underlying polymer configurations can eas- 
ily extend to distances greater than 2Rg, the POA, which 
implicitly assumes that the interaction of an effective par- 
ticle with one wall is not directly affected by the presence 
of a second wall, begins to break down for strong confine- 
ment. Interestingly, the trend shown in Fig. G4 for the 
CM effective particle representation of the Gaussian coils 
mirrors the trend shown in Fig. 11^ for interacting SAW 
polymers, suggesting that the differences for the latter 
also arise from the breakdown of the POA approxima- 
tion. 

We note that end-points or mid-points could also be 
used to construct an effective particle picture of interact- 
ing polymer solutions. For example, the mid-point rep- 
resentation would be very similar to the two-arm limit 
of a star-polymer, for which a number of results have 
been recently derived [p7| , p9|j6^ ] . There are differences 
with the CM representation; for example, the mid-point 
equivalent of Eq. (0) would scale as |^,0 : 



V2ir) 



9 \R 



(A9) 



at short distances. Here i? is a lengthscale proportional 
to the polymer size Rg. In contrast to the CM case, 
this interaction diverges at full overlap so that one would 
expect some qualitative differences in the behavior of 
the underlying soft particle fluids. However, the two 
approaches should, in principle, produce similar results 
for the thermodynamic properties of polymer solutions. 



The relative merits of using mid-points v.s. the CM to 
describe polymer solutions are currently under investiga- 
tion. 

Finally we note that several results for ideal polymers, 
obtained here by using end-point or mid-point densities, 
can be also obtained by using the exact monomer den- 
sity profiles near one wall U^ together with the DOA 
approximation 173]. 
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